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I. INTRODUCTION 



The analysis of the flavour structure of the neutrino mass matrix has been a subject of 
intensive study recently |L|]. A non-trivial neutrino mass matrix can solve, through the os- 
cillation solutions, the atmospheric and solar neutrino anomalies, which have been observed 
by different experiments: to the very recent results of the Super-Kamiokande collaboration 
0, have to be added those of other atmospheric neutrino experiments (1MB ||, Soudan 
||], Kamiokande 0) and solar neutrino experiments (Homestake ||, Gallex J7|, SAGE ||, 



Kamiokande [|, Super-Kamiokande |g, MACRO JTTJ and LSND pj 

The oscillation explanation of the solar neutrino problem, the atmospheric neutrino 
anomaly and the LSND results suggest three very different values of neutrino mass squared 
differences, namely Am s 2 un < Ara 2 atB < Am£ SND , with Am^ < 10" 4 eV 2 , Am 2 tm G 
[10~ 3 , 10~ 2 ] eV 2 and AmL SND G [0.3, 1.0] eV 2 . The evidence for — > v e observed by 
LSND has not been confirmed or excluded by the KARMEN [^] experiment. MiniBooNE 
at FNAL fl4| or MINOS long-baseline experiments |15] could provide the answer. Excluding 
LSND results, the oscillation explanation of the solar neutrino problem and the atmospheric 
neutrino anomaly requires two mixing angles and two values of neutrino mass squared dif- 
ferences with a strong hierarchy, namely Am 2 un <C Am 2 tm . 

In a previous analysis |16[ with a general (3 x 3) symmetric mass matrix assuming 
oscillations as an explanation of the neutrino deficits, we had used as inputs neutrino data 
constraints in order to obtain the allowed range of variation of the mass matrix elements. 
These inputs are Am 2 tm , sin 2 2# atm from SuperK measurements, Am 2 un , sin 2 26* sun from solar 
experiments and finally the CHOOZ constraint sin 2 6'cHooz- We also discussed the 
spectrum after diagonalization and considered two interesting quantities: the sum of the 
neutrino masses, which is a relevant quantity for obtaining a neutrino component of hot dark 
matter (HDM), and the effective mass constraint, which is relevant for neutrinoless double 
beta decay (this is applicable when the neutrinos are Majorana and massive particles). 

To generate a Majorana mass for the neutrinos in a given model lepton number must 
be violated. We consider the Minimal Supersymmetric Standard Model (MSSM) with R- 
parity \TE\ violation (Rp = (_i) i + 3 - B + 2S ) where L, B, S are the lepton and baryon number 
and the spin of the particle, respectively), left-handed neutrinos obtain a Majorana mass, 
at tree level, through mixing with the neutralinos, and through loop diagrams that violate 
lepton number (in two units). The simultaneous presence of baryon and lepton number 
violating couplings is not acceptable, because of the long lifetime of the proton. As we are 
interested in a model that provides lepton-number violation, we simply choose all baryon 
number violating couplings to be exactly zero. 



In ref. |16[ we presented as an application of our results the case of the Minimal Super- 
symmetric Standard Model (MSSM), with R-parity violation at the one- loop order allowing 
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the presence of both bilinear and trilinear Rp-violating couplings. We included the effect 
of the bilinear couplings at tree level and of trilinear couplings at one loop. The numerical 
bounds were given assuming all trilinear \j k = A, as well as A^ fe = A' to be equal. Here we 
will relax this assumption and apply our general results to other limiting (albeit general) 
cases in this model. The aim is to provide strong constraints on the couplings with an ex- 
plicit reference to the generation indices from neutrino data. We will allow the simultaneous 
presence of various lepton number violating couplings. 

The paper is organized as follows. In section II we describe the inputs used to constrain 
the R-parity violating parameters, while in section III we introduce our notation for the 
MSSM with Rp-violation. We also give in section III the tree-level neutrino (3 x 3) mass 
matrix and the loop corrections to each matrix element. We parametrize the mass and 
mixing matrices. In section IV, we give results for five separate subcases that correspond to 
some limiting cases of the MSSM without R-parity, and we derive bounds on the bilinear 
and trilinear Rp-violating couplings for each case. 



II. GENERAL INPUTS 

Global fits of neutrino data have shown that neutrino oscillations among three flavours 



are sufficient to accommodate the solar and atmospheric data [19]. Three different types 
of solar experiments are sensitive to different solar neutrino energy ranges; consequently, 
three different ranges of solutions for the solar data exist, which correspond to the vacuum 
oscillation solution, MSW with a large mixing angle (MSW-LMA) and MSW with a small 
mixing angle (MSW-SMA). The required neutrino mass squared differences and mixing 
angles are shown in table I. We use these inputs together with the results from the CHOOZ 



experiment |Ty| to constrain the elements of our neutrino mass matrix in the flavour basis. 
An additional sixth input is needed to solve the general (3 x 3) real and symmetric mass 
matrix. As our sixth input we could take the direct upper bound on the effective mass from 
neutrinoless beta decay (/3/?)o„, or the upper bound on the effective neutrino mass coming 
from the 3 H /3-decay spectrum, or the astrophysical bound on the magnetic moment of the 
neutrino bounds, or the global cosmological upper bound on the sum of the neutrino masses 



Let us briefly discuss each of these possible inputs separately. The direct upper bound 
from the measurements of the high energy part of the 3 H /3-decay spectrum given as upper 



limits on the electron neutrino mass obtained in the TROITSK [21] and MAINZ [22] exper- 
iments are m v < 2.5 eV and m u < 2.8 eV, respectively. However, these experiments suffer 
from some ambiguities referred to as "the negative mass squared problem", which is still 
not completely understood and we will therefore not use these bounds. 
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The upper bound on the sum of the neutrino masses coming from astrophysical and 
cosmological considerations, such as the one from hot dark matter (HDM), suggests that 



m n = mi + m 2 + m 3 < few eV, 

i 

where rrii,i = 1,2,3, are the masses of the mass eigenstates z/j, which constitute active 
flavour neutrinos. A recent analysis, based on the observation of distant objects favouring 
a non-zero cosmological constant instead of HDM [p3|l , shows that the bound of a few eV 
for the sum of the neutrino masses is no longer required, since a HDM component is not a 
necessary ingredient in this case. 

In the R-parity violating model that we are considering we can have transition magnetic 
moments that change lepton number. These magnetic moments are proportional to the 
matrix elements; we prefer to use as an input a direct bound on a matrix element as such 
arising from neutrinoless double beta decay. 

The additional input we will take is the best limit on the effective mass m e g appearing 
in (/5/?)o^, as it directly constrains the matrix element mu of the neutrino mass matrix. It 
is defined by 

Kff| = \ J2 m »i U ei\ <Y, m »k\ U ekl (!) 
i k 

which has been derived in the Heidelberg-Moscow 76 Ge experiment |24| (see also references 
I2R 



|m eff | < (0.2 - 0.6) eV (at 90%CL). (2) 

The CP phases that might appear in the mixing matrix elements U e i are not relevant to our 
analysis, as we only use the second r.h.s. term in eq. ([I]). 

For three generations, the flavour states V\ are expressed in terms of the mass eigenstates 
Vi using the (3 x 3) mixing matrix U 

i/, = £l/Ji/ 4 . (3) 
i=i 

Using the Chau and Keung parametrization of a 3 x 3 rotation matrix ||27|, U is given by 



U 



Cl2 C13 S\2 C13 S\ 3 e \ 

-Sl2 C 23 - C12 S 2 3 Si 3 e %S C 12 C 23 - Sl2 S 23 Si 3 e lS S 23 C 13 
V S12 S 23 - C12 c 23 si 3 e lS -c 12 s 23 - s 12 c 23 Si 3 e tS C 23 C13 / 



diag{e iai ,e iQ2 ,l}, (4) 



where Cjj = cos(^j) and Sy = sin(^j), 5 is the Dirac CP phase and a\ y2 are the Majorana 
ones (we have two additional CP phases in the case of Majorana particles). 
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The survival probability P„ e _M/ e , relevant to the case of solar fluxes, depends only on the 
first row of the mixing matrix in eq. @, i.e. on |£/ e j| 2 , with i = 1,2,3. In the atmospheric 
case, the oscillation probability depends on the last column of i.e. on jc/^ 2 , with 
£ = e, ix, t. The other elements of the matrix are not constrained by any direct experimental 
observation. With this parametrization, we directly obtain the mixing angles sin 2 26\2 and 
sin 2 2#23, and also the relevant CHOOZ parameter sin #13. 

We impose the following hierarchy on the mass eigenvalues and denote them by m^, i — 
1, 2, 3, such that m\ < m 2 < 7713. As the neutrino oscillation scenario cannot fix the absolute 
mass scale nor distinguish whether the smallest mass splitting is between the two lightest 
mass eigenstates or the two heaviest ones, we present in table II the possible neutrino spectra 
and indicate which is the corresponding mass squared difference. 

To summarize, we will use the following six inputs from neutrino data: 

• the limit from the effective mass appearing in neutrinoless double beta decay: m e g, 

• Am 2 tm and sin 2 20 atm , 

• Am 2 un and sin 2 26» sun , 

• sin 2 2<9 C Hooz- 



Experiment 


Am 2 (eV 2 ) 


sin 2 29 


Atmospheric 


(2 - 5) x 1(T 3 


0.88 - 1 


Solar 






MSW-LMA 


(3 - 30) x 10~ 5 


0.6 - 1 


MSW-SMA 


(0.4-1) x 10~ 5 


10" 3 - 10~ 2 


Vacuum 


(0.5 - 8) x 10~ 10 


0.5 - 1 


CHOOZ 


> 3 x 10~ 3 


< 0.22 



TABLE I. MSW-LMA, MSW-SMA and Vacuum stand for MSW large mixing angle, small 
mixing and vacuum oscillation solutions, respectively. 
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Spectrum 


Solar 


Atmospheric 


Hierarchy 


Am 2 12 


Amf 3 


Degenerate 


Am^ or Amf 2 


Amf 3 


Pseudo-Dirac 


Am^ or Awij 2 


Amf 3 



TABLE II. Different possible regimes and corresponding mass squared difference. 



III. MSSM WITH R P - VIOLATION 



The most general renormalizable superpotential for the supersymmetric Standard Model 
with lepton-number violation is 

W = e ah {n a L a a H h u + \ a p k L%L^ + h? k QtH»Ug + \' aik LiQfD^, (5) 

where the (i,j,k) are flavour indices, (a, b) are SU{2) indices, and the (a, (3) are flavour 
indices running from to 3. The L a are the doublet superfields with hypercharge Y = — 1. 
Note that the A couplings are antisymmetric in the first two indices. The usual R-parity- 
preserving Lagrangian is obtained when only Aoifc = hf k ,X' i0k = hf k are non-zero, and 
we can identify L = Hj. In the model of eq. (|5|) we have 9 additional A couplings and 
27 new A' couplings with respect to the R-parity-conserving case. Note that thanks to the 
additional degrees of freedom, we can rotate in the flavour space of the "down-type" scalar 
fields to set the vacuum expectation values of the sneutrinos to be zero Q. Henceforth, we 
will work consistently in this basis and the bounds we will derive are valid for this basis. 
This model has been extensively analysed in the literature p8||- [[31]. References |30|j31 



showed in a basis-invariant way that neutrino masses are always generated in these models, 
even when universality of the soft SUSY-breaking terms is assumed at some high scale. 
Previous studies have also tried to constrain the different Rp-violating couplings that appear 
in the MSSM Lagrangian, considering only the effect of bilinear terms [B^] or only of trilinear 



couplings |[33|| , or of both [|34|,|T6|| , from solar and atmospheric neutrino data. Both tree- 



level and one-loop effects have been considered. A recent study [^5] has constrained these 



couplings using rare decays. The results are in agreement with the constraints given in ref. 



It is well known that the tree-level expression of the neutrino mass matrix elements 
obtained through the mixing with neutralinos is 



This can be done order by order in the loop expansion when one appropriately defines the mass 
matrices of the Higgs sector. 
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slept ons 



leptons 



squarks 



quarks 



FIG. 1. One-loop Feynman diagrams contributing to the neutrino masses. 



M 



tree 



2 (Mi + tan 2 9 W M 2 ) 2 
92 4delM Wl 



9l 



v 2 cos 2 j3 
M 3 

susy 



(6) 



where i, j are flavour indices, Mi, M2 are the gaugino masses and det M is the determinant 
of the R-parity conserving neutralino mass matrix. For the second term of the right-hand 
side we have taken all R-parity conserving masses to be of order ~ M susy and as we are 



working in the basis where the slepton vacuum expectation value is zero, v 1 = Vd 
We can rewrite the tree-level mass matrix as 



M 



tree 



c 



( H 2 e fi,j,fl e {l e /J T \ 



v cos (3. 



(7) 



The above mass matrix has only one non-zero eigenvalue. In order for the other neutrinos 
to obtain a mass, we must consider loop corrections. There are many types of loop diagrams 
that can contribute to the neutrino mass in this model [29 , 55 , ^7 . Here we will focus on the 
contribution at one loop from fermion-sfermion diagrams involving trilinear couplings, see 
fig. [[} See, for example, ref. [pl| for more details regarding the computation of the neutrino 
mass matrix. 

The one-loop mass contribution from slepton-lepton and squark-quark loops to each 
element of the mass matrix is given by 



qm ) loop 



16tt 2 



X 



J( X q 



Ml 



M h k, P 



E\i \> ... 
"qkp^mpk 



(8) 
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where 

where X is the trilinear term that appears in the off-diagonal matrix element of the sfermions 
mass matrix. We consider here that we are in the down-quark mass eigenstate basis, and 
that the A, A' couplings have been redefined in terms of the couplings appearing in the 
superpotential and of the corresponding Cabibbo-Kobayashi-Maskawa matrix elements; we 
also drop the flavour dependence from 2 /(x^) and 2 f(x^), and consider them to 
be universal in the slepton and in the squark sector, respectively. Thus the matrix elements 
at one loop are given by the sum of eqs. ([F]) and (|8|). The analysis we perform can directly 
constrain terms of the form Cfiifij, K[ J J2k, P Kkp^jpk, K'2 I2k, P KkpKpk> wriere 



X f(x e ) 
16tt 2 M? 



K i = T^^HPT ( m i m j) > 



^ 3 i§^r^- < 10 > 

Note that the quantities K{ 3 2 depend only on R-parity conserving parameters. As an ap- 
plication, we focus in the following on a case where we assume that the A- jfc are of the 
same order, which allows us to neglect the extra mixed A' 123 132 223 232 due to the hierarchy 
m s mb <C ml (see eq. (§)). We also assume separately that all Xijk are of the same order 0. 
Thus, using the mass hierarchy m e> d <C rn^ tS <C ra T) & in (m gm )i oop , simplifies the expressions 
for the matrix elements, to obtain 



M 



loop 



( ^1^133 + -^Ailg -^iAi3 3 A 2 33 + -^2A' 133 A 233 K 2 A' 133 A 333 \ 

-^lAi33A 23 3 + ^2Ai 33 A 2 33 i^lA 233 + K 2 X 2 ^ ^2A 233 A 333 



i/2 



(11) 



\ -^2A' 13 3A 333 ^2^233 A333 ^2^333 / 

where the coefficients are given by: 

X f( Xt ) ( 2 \ 



Thus, we take our one-loop mass matrix to be 



2 The same procedure as we apply here could be performed if we considered that only one of the 
R-parity violating trilinear couplings was dominant. 
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M v = M^ 00 + M 1 ° op . (13) 

At this order, the mass matrix has in general three non-zero eigenmasses. The 8 R-parity 
violating parameters that we would like to constrain are: 

• ^133) ^233 

• ^133 ; ^233 an d A 333 

• fi e , fi^i and fi T . 

Recall that we will use only 6 inputs from solar (2), atmospheric (2), CHOOZ (1) constraints 
and the bound on the effective neutrino mass appearing in neutrinoless double beta decay. 
In our previous analysis fl6| , denoted case in table III, we considered a toy model 



leading to three non-zero mass eigenstates. In this toy model we assumed that all trilinear 
couplings Xijk were equal and that all A^, fe couplings were equal, i.e. = A and X'^ k = A'. 
This leads to five unknowns f^ e ,n,r > A and A', and the system was solved using the five inputs 
from the solar, atmospheric and CHOOZ constraints. 

In this work, we will relax this assumption to more general cases. We consider five 
different cases that generate three non-zero eigenmasses. These different subcases of the 
generic matrix of eq. (|T3|) , arise when we apply certain conditions on some of the L-number 
violating parameters. We do this so as to reduce the number of unknowns to the number of 
constraints. The subcases we consider are summarized in table [| 

A first case corresponds to the situation where all the A' couplings are switched off, 
which happens in the limit in which the squarks decouple. Case 2 occurs when all the A 
couplings are switched off, i.e. in the limit in which the sleptons decouple. In case 3, the 
bilinear couplings are such that [i e ~ [i^ <C \i T and thus only one bilinear coupling (/i r ) is 
relevant. In cases 4 (5), the trilinear A (A') are of the same order, i.e. A133 ~ A233 = A 
(A' 133 ~ A'233 ~ A333 = A') and the bilinears are such that fi e ~ /i M 7^ /j, t for case 4 and 
A*e 7^ /V 7^ f° r case 5- We do not claim to have covered all possibilities, but certainly 
some of the most straightforward ones. 
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Cases 


Kjk 


Kjk 


Hi 





Kjk = A 


Kjk = a' 


Me> M/jj Mr 


1 


Al33, A233 





Me) M/^ Mr 


2 





A 133' A 233' A 333 


Mej M/ij Mr 


3 


<Vl33i ^233 


A 133' A 233' A 333 


Mr 


4 


Kjk = A 


A 133' A 233' A 333 


Mi">Me — M/i 


5 


^133) ^233 


Kjk = k 


Me; M/i) Mr 



TABLE III. Different cases of Rp-violating couplings contributing to the neutrino mass matrix. 



IV. RESULTS 



In our general scan of parameter space, we allow tree-level contributions to either domi- 
nate over the loop corrections, or to be on the same order as these, or to be much smaller 
than the loop terms. The analysis for subcase 0, presented in fl6|], gave bounds on the R p 



couplings from combinations of constraints from atmospheric and CHOOZ data, together 
with one of the possible solar neutrino solutions. For our five new subcases we relax the 
assumptions of subcase to present stringent constraints on bilinear and trilinear lepton- 
number violating couplings, with specific generation indices in the basis where the sneutrino 
vacuum expectation value is zero. 

The bounds on the couplings are presented in tables for the various combination of 
constraints from neutrinoless double beta decay, atmospheric, CHOOZ together with vacuum 
or MSW-SMA or MSW-LMA solutions. A common assumption used to place bounds in this 
model is to take all Rp-conserving mass parameters to be of the same order, M susy . For the 
particular case where M susy = 100 GeV, f(x) — > 1 and tan/3 = 2, we have 

K x ~ 1.8 x 10~ 4 GeV, K 2 ~ 4.7 x 10~ 3 GeV, C ~ 5.3 x lO^GeV -1 . 

We use these values to obtain the numerical results presented for cases 1-5 in tables IV- VIII, 
respectively. Using eqs. (|) and ( |T2"D modified bounds can be obtained for other values of 
the R-parity conserving parameters. 

The analysis for all subcases shows that there are many regions of parameter space that 
can simultaneously accommodate the MSW solution (large or small mixing angle), SuperK 
and CHOOZ constraints. Subcases 2-4 also provide solutions for the combined constraint 
of atmospheric data, CHOOZ and the vacuum oscillation solution, while subcases 1 and 5 
do not. The results for the latter two subcases are similar to those obtained for subcase 
treated in [0. The number of solutions defining the allowed region in the (Am 2 , sin 2 20) 
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Couplings 


71 THTTT TATA 

MSW-LMA 


TV THTTT fl TV /T A 

MSW-SMA 


Vacuum 


1-^133 1 


6.8 x 1CT 4 


1.3 x 1CT 4 




AoQQ 
| xv ZOO | 


5.6 x 1CT 4 


1.9 x 1CT 4 




W\ (GeV) 


6.1 x 1(T 5 


2.o x i(r 6 




W\ (GeV) 


1.2 x 10~ 4 


1.2 x 1(T 4 




iMrl (GeV) 


1.2 x 1(T 4 


1.2 x 10 -4 





TABLE IV. Bounds on the couplings for subcase 1 that satisfy MSW-LMA or MSW-SMA, 



SuperK, CHOOZ and neutrinoless double beta decay simultaneously. 



Couplings 


MSW-LMA 


MSW-SMA 


Vacuum 


1 \' 1 

A 133l 


1.5 x 10~ 4 


4.1 x 1(T 5 


5.4 x 1(T 5 


1 X' 1 
\ A 233\ 


1.5 x 10" 4 


1.5 x 10" 4 


9.5 x 1(T 5 


1 \' 

\ A 333l 


1.5 x 1(T 4 


1.5 x 1(T 4 


1.5 x 1(T 5 


lMe| (GeV) 


1.2 x 10~ 4 


4.1 x 1(T 5 


5.1 x 1(T 5 


M (GeV) 


1.2 x 10" 4 


1.2 x 1(T 4 


8.1 x 1(T 5 


\fir\ (GeV) 


1.2 x 10" 4 


1.2 x 1(T 4 


1.2 x 10~ 4 



TABLE V. Bounds on the couplings for subcase 2 that satisfy MSW-LMA or MSW-SMA or 
vacuum, SuperK, CHOOZ and the neutrinoless beta decay constraints simultaneously. 



plane is very small for the vacuum solution, larger in the case of MSW-SMA, and it is still 
larger in the case of MSW-LMA. This is in agreement with the most recent results from the 
Super-Kamiokande collaboration, presented at the SUSY2K conference [ftBfl , stating that 
the solar neutrinos favour the MSW-LMA solution. There is less room for MSW-SMA 
and much less still for vacuum solutions. We can see that our bounds are consistent with 



bounds derived in the literature (see p9j-44~T| for three-neutrino case and |^2[ for four-neutrino 
case). We emphasize that some of our subcases present more stringent bounds, and that 
the combination of constraints with the vacuum solution requires the smallest values for 
R-parity violating couplings. 

To summarize, we have obtained bounds on the bilinear and trilinear R-parity violating 
couplings with explicit reference to the leptonic indices. The bounds have been obtained 
from neutrino data, which constrain the neutrino mass matrix that can be constructed in 
the MSSM with R-parity violation for three generations of neutrinos. We considered the 
tree-level contribution and the one-loop contribution from fermion-sfermion diagrams with 
trilinear couplings to the neutrino mass matrix. 
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Couplings 


MSW-LMA 


MSW-SMA 


Vacuum 


1 \ f 1 

l A 133l 


1.2 x 1CT 4 


6.8 x 1CT 5 


1.4 x 1CT 5 


1 \ t 1 

l A 233l 


1.5 x 10 -4 


1.5 x io~ 4 


1.4 x 1CT 5 


1 V 1 
1^333 1 


1.5 x 10~ 4 


1.5 x 10~ 4 


1.5 x 10~ 4 


^133 


6.2 x 1(T 4 


3.1 x i(r 4 


1.1 x 1(T 7 


^233 


6.2 x 10~ 4 


3.1 x i(r 4 


1.1 x 1(T 7 


lA*r| (GeV) 


1.3 x 10 -4 


1.3 x 1(T 4 


1.2 x 1(T 5 



TABLE VI. Bounds on the couplings for subcase 3 that satisfy MSW-LMA or MSW-SMA or 
vacuum, SuperK, CHOOZ and the neutrinoless beta decay constraints simultaneously. 



Couplings 


MSW-LMA 


MSW-SMA 


Vacuum 




1.1 x 1(T 4 


3.4 x 1CT 5 


5.4 x 1(T 5 


^233 


1.4 x 1(T 4 


1.4 x i(r 4 


5.4 x 1(T 5 


l A 333l 


1.5 x 1(T 4 


1.4 x i(r 4 


1.5 x 1(T 4 


|A| 


4.4 x 10~ 4 


1.9 x 10" 4 


1.1 x 1(T 7 


\n e ~ (GeV) 


8.1 x 1(T 5 


2.3 x 1(T 5 


4.6 x 1(T 5 


\Hr\ (GeV) 


1.2 x 1(T 4 


1.3 x 1(T 4 


1.3 x 1(T 4 



TABLE VII. Bounds on the couplings for subcase 4 that satisfy MSW-LMA or MSW-SMA or 
vacuum, SuperK, CHOOZ and the neutrinoless beta decay constraints simultaneously. 
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Couplings 


TV imi7 T TV /T A 

MSW-LMA 


MSW-SMA 


Vacuum 


1-^133 1 


6.2 x 10~ 4 


1.3 x 10~ 4 




I A233 1 


5.6 x 1CT 4 


3.1 x 1CT 4 




|A'| 


6.8 x 10~ 5 


4.1 x 10~ 5 




W\ (GeV) 


8.1 x 1(T 5 


3.4 x 1(T 4 




M (GeV) 


1.3 x 1(T 4 


1.3 x 10 -4 




Kl (GeV) 


1.3 x 10~ 4 


1.3 x 1(T 4 





TABLE VIII. Bounds on the couplings for subcase 5 that satisfy MSW-LMA or MSW-SMA 
or vacuum, SuperK, CHOOZ and the neutrinoless beta decay constraints simultaneously. 
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